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ABSTRACT 

The quasar SDSS J133401.39+331534.3 at z = 2.426 is found to be a two- 
image gravitationally lensed quasar with the image separation of 0'.'833. The 



rj I object is first identified as a lensed quasar candidate in the Sloan Digital Sky 



Survey Quasar Lens Search, and then confirmed as a lensed system from follow- 



On' up observations at the Subaru and University of Hawaii 2.2-meter telescopes. We 

O ■ estimate the redshift of the lensing galaxy to be 0.557 based on absorption lines in 



the quasar spectra as well as the color of the galaxy. In particular, we observe the 
system with the Subaru Telescope A0188 adaptive optics with laser guide star. 
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in order to derive accurate astrometry, which well demonstrates the usefulness 
of the laser guide star adaptive optics imaging for studying strong lens systems. 
Our mass modeling with improved astrometry implies that a nearby bright galaxy 
~ 4" apart from the lensing galaxy is likely to affect the lens potential. 



Subject headings: gravitational lensing 
— instrumentation: adaptive optics 



quasars: individual (SDSS J133401. 39+331534.3) 



Introduction 



Since the discovery of the first gravitationally lensed quasar (j Walsh et al.lll979[ ). it has 
convincingly been demonstrated that lensed quasars provide insights into various research 
fields in astrophysics and cosmology, as well as being a unique tool for studying the dark 
universe. For instance, one can directly measure the structure and substructure of lensing 
objects, including the distri bution of dark matter, through detailed observations and anal- 



ysis o f lensed systems (e.g., iTreu fc KoopmansI I2004J : iRusin fc Kochanekl l2005t IChiba et al. 



20051 ). Magnification by gravitationa l lensing greatly enhances our ability to study quasar 



host galaxies at high-redshifts (e.g., iPeng et al.l l2006l ). Moreover, the statistics of lensed 



quasars, as well as measurements of time del ays between lensed images, provide independent 
constraints on co smological parameters (e.g. 
Suvu et al.ll201ol ). 



Chae et al.ll2002l : IOgurill2007l : lOeuri et al.ll2008 



High resolution imaging such as that provided by Adaptive Optics (AO) is crucial for 
accurately characterizing lensed quasar systems. This is not only the case for the small- 
separation systems (~ 1"), where accurate astrometry (prone to the effect of atmospheric 
seein g) is necessary for constraining mass models of the lensing galaxies (e.g., ISluse et al. 
20081 ). but also for the larger separation systems, where adaptive optic s imaging can identif y 



lens substructures and faint extended lensed features such as arcs (e.g.. lMcKean et al.ll2007l ). 



The Sloan Digital Sky Survey Quasar Lens Search (SQLS; lOguri et al.l l2006l 12008 



Inada et al.ll2008l l2010l ) has been one of the most successful strong lens surveys conducted to 
date. SQLS relies on the large homogeneo us sample of spec troscopically-confirmed quasars 
from the Sloan Digital Sky Survey (SDSS; lYork et al.l 120001 ). By identifying lensed quasar 
candidates using the color and morphological information from the SDSS imaging data, 
SQLS has discovered more than 40 new lensed quasars so far. Together with previously 
known lenses which were re-observed by the SDSS, the SQLS sample consists of ~ 60 lensed 
quasars, constituting roughly half of all lensed quasars discovered to date. 



However, a disadvantage of SQLS, like other ground-based optical strong lens surveys. 
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is its poor capability of identifying subarcsecond (image separation 9 < 1") lensed sys- 
tems. This is because typical seeing sizes of ground-based optical imaging observations are 
~ 1", which makes it difficult to resolve multiple components of subarcsecond systems. In- 
deed, discoveries and c onfirmations of suba rcsecond lensed quasars have been made mostly 
by radio surveys (e .g., iBrowne et al.l 120031 ) or by using the Hubble Space Telescope (e.g.. 



Reimers et al. 



2006 



20021). with few exceptions from ground-based surveys (e.g.. ICastander et al. 



Blackburne et al.ll2008h . 



In this paper, we report the discovery of the first new subarcsecond (two-image) lensed 
quasar discovered in the course of SQLS, SDSS J133401. 39-1-331534. 3 (hereafter SDSS J1334-I-3315). 
After briefiy describing the lens candidate selection from the SDSS data (§2]), we present our 
imaging and spectroscopic follow-up observations (§3]). We further observe this system with 
the Subaru Telescope laser guide star adaptive optics (LGS-I-A0188), for robust detection 
as well as characterization of the lensing galaxy (^U ^. The accurate astrometry obtained 
with the LGS-I-A0188 is then used for gravitational lens mass modeling (^. We summarize 
our results in ^ We assume the concordance cosmology with Hq = 70 km~^ s~^ Mpc~^, 
Qm = 0.27 and ^a = 0.73 throughout this paper. 



2. Candidate Selection from the Sloan Digital Sky Survey 



The SDSS (SDSS-I and SDSS-II Legacy Survey; lYork et al.l 120001 ) is a combination of 
imaging and spectroscopic surveys to map 10,000 square degrees of the sky, cen tered at the 



North Galactic Cap. The survey uses a dedicated wide-field 2.5-meter telescope (JGunn et al. 



20061 ) at the Apac he Point Observatory in New Mexico, USA. Images ta ken in five broad 
band filters (ttori^^. lFukugita et al.lll996l : iGunn et al.lll998l : iDoi et al.ll2010f ) are re duced with 



an au tomated pipeline, achieving astrometric accuracy better than about O'.'l (jPier et al. 



20031 ) . and photometric zeropoint accurac y better than about 0.01 magriitude over the en- 



tire survey area, in the q, r, and i bands (IHogg et al.l l200ll : ISmith et al.l l2002l : llvezic et al. 



2004J : iTucker et al.l 120061 : iPadmanabhan et al.ll2008l ). The SDSS also conducts spectroscopic 



observations with a multi-fiber sp ectrograph covering 3800 A to 9200 A, with a resolution 



of i? ~ 1800 fBlanton et al. 



2OO2I : lAbazajian et al. 



2003 



2004 



2003 ) . The data are made publicly available (IStoughton et al. 



20051 . I2OO9I : lAdelman-McCarthv et al.ll2006 



2007 



20081 ). 



The SQLS (lOguri et al.ll2006l . I2OO8I : llnada et a 



20081. 



20101) identifi es lensed quasar 



can- 



didates among spectroscopic quasars in the SDSS (JRichards et al.ll2002l ). using two selection 
algorithms. One is the morphological selection, which selects quasars with possible extended 
morphology as small-separation lensed candidates. The other is the color selection, which 
searches for nearby stellar objects whose colors are similar to the spectroscopically identified 
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quasars. SDSS J1334+3315 is a lensed quasar candidate a.t z = 2.426 selected by the mor- 
phological selection (see Figured] for a finding chart). Although the morphological selection 
is designed to be nearly complete at ^ > 1", it can identify some sub arcsecond leiised sy stems 



as well, as explicitly shown in simulations of lensed quasar images ( lOguri et al.ll2006[ ). The 
point spread function (PSF) magnitudes of the unresolved SDSS J1334+3315 system, mea- 
sured by the SDSS, are u = 19.32 ±0.03, g = 18.76 ±0.01, r = 18.66 ±0.02, i = 18.80 ±0.02, 
and z = 18.71 ± 0.03, where the Galactic extinction is not corrected. 



3. Follow-up Spectroscopy and imaging 
3.1. Spectroscopy 

We conducted spectroscopic follow-up o bservations of SDSS J1 334+3315 with the Faint 



Object Camera and Spectrograph (FQCAS: iKashikawa et al.ll2002l ) at the Subaru 8.2-meter 



telescope on January 31, 2009. We adopted the 300B grism and the L600 filter to cover the 
wavelength from 3700 A to 6000 A. We obtained data with an exposure of 600 s using a 1" 
slit, with 2x2 on-chip binning. The resulting spatial resolution is 0'.'208 pixel"^, and the 
spectral resolution is R ^ 400. The observations were necessary in order to extract individual 
spectra of the two stellar components (the two quasar images), which are unresolved in the 
original SDSS spectroscopy. The one-dimensional spectra of the stellar components were 
extracted using standard IRAF tasks. We were unable to extract the spectrum of a lensing 
galaxy (see §3.21) . due to its faint nature and the subarcsecond separation of the quasar 
components. Although the two stellar components are resolved reason ably well in the two- 



dimensional spectrum, we adopted a deblending procedure developed by lPindor et al.l (120061 ) 
in extracting one- dimensional spectra, in order to minimize the cross-contamination between 
the two components. 

The spectra shown in Figure [2] indicate that both stellar components have quasar broad 
emission lines (Lya, Si IV, C IV) at the same wavelengths, corresponding to a redshift of 
2.426 (measured also by the initial SDSS spectroscopy). The emission lines have the same 
overall shape, supporting the gravitational lensing hypothesis. Several absorption lines are 
also observed in both spectra, at the same wavelengths. We identify some of these as the 
2599 A Fe II and the 2796 A, 2803 A Mg II doublet, from which we derive the redshift of the 
absorber to be 2; = 0.557. The absorption is seen in both components, with the absorption in 
component A being slightly stronger than that in component B. We postulate this absorber 
to be the lensing galaxy, as the redshift is consistent with the photometric redshift estimate 
of the lensing galaxy, as we will discuss later (^. In addition, we identify a strong C IV 
absorber at z = 2.172 in both stellar components. 
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The ratio of the fluxes of the two stellar components is quite constant along the plotted 
range (B/A = 0.7 - 0.8). There is therefore no sign of differential dust reddening which, 
should it occur on the line of sight to one of the quasars, would affect the flux ratio towards 
shorter wavelengths. Because the flux ratio is quite constant both across the continuum and 
the broad line regions, there is also no evidence of microlensing effects. Microlensing due 
to stars in the lensing galaxy would effect the continuum emission more strongly than the 
broad line emission, since the physical size of the regi on associated wit h the former is much 
smaller than the one associated with the latter (e.g., ISluse et al.l 120081 ). The change in flux 
ratio below ~ 4000 A in the observer's frame is thought to be due to errors arising from the 
low detector response in that region. 



3.2. Optical and Near-Infrared Imaging 

Optical (J and z) follow-up images were taken with the Tektronix 2048x2048 CCD 
camera (Tek2k) at the University of Hawaii 2.2-meter (UH88) telescope on 2009 April 15. 
The pixel scale is 0'.'219 pixel"^. The total exposure time was 800 s for each filter. The seeing 
was ~ 1" — 1'.'2. The images were reduced using standard IRAFj^ tasks. The photometric 
zero-points were derived by comparing magnitudes of nearby stars located in the same frame 
to their SDSS magnitudes, with the resulting zero-point accuracy of ~ 0.05 mag. For the / 
band, the magnitude conversion has been performed using the equations of Lupton 2005o 

In addition, we obtained near-inf rared (JHK,^) images with the Multi- O bject InfraRed 
Camera and Spectro graph (MOIRC S; llchikawa et al.ll2006t ISuzuki et al.ll2008l ) at the Subaru 
8.2-meter telescope ( llye et al.ll2004j ) on 2010 April 2, under an excellent seeing size of ~ 0'.'4. 
The pixel scale is 0'.'117 pixel"^. The total exposure time was 600 s for J and Ks bands, and 
480 s for H band. We use the MCSRED software package (I. Tanaka et al., in preparation) 
for the data reductio n. We derived the magnitude zero-points using the standard star FS33 
( JLeggett et al.l 120061 ). leading to a zero-point accuracy estimated at ~ 0.05 mag. 



Results of the follow-up imaging are summarized in Figure El In addition to the can- 
didate lensed quasar system, there are two red galaxies (Gl and G2) within 5". In order 
to check the presence of the lensing galaxy in between the stellar components, as well as 
to extract relative astrometry and photometry from the data, we analyze the images using 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 



http : //www . sdss . org/drS/algorithms/sdssUBVRITrcinsf orm . html 
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the public software GALFIT (jPeng et al.ll2002l ). We fit all components simultaneously with 
nearby stars as PSF templates and galaxies modeled by the Sersic profile. For the bright 
galaxy Gl, the preliminary GALFIT fitting produces the best-fit Sersic index oi n = 0.87, 
very close to the canonical value n = 1 for disk-like galaxies. The faint galaxy G2 appears to 
be quite red, and is also not detected in the optical images. Due to its faintness, we cannot 
extract an accurate Sersic index. 

Table [U provides a summary of the GALFIT modeling results. The separation between 
the stellar components (the two quasar images), derived from the Kg-hand image, is 0'.'83 ± 
O'.'Ol. We detect an extended object (galaxy G) in between the two stellar components (the 
two quasar images), which can be interpreted as a lensing galaxy. We find that the lensing 
galaxy is hardly detected in the UH88 optical images, while it is significantly detected in the 
Subaru/MOIRCS near-infrared images (see also Figure [3]). The difference can be attributed 
to the much better angular resolution and the deeper depth of the MOIRCS images. The 
position of the lensing galaxy G is difficult to determine accurately, as the best fit varies 
between different bands. However, by adopting the position estimated in the i^s-band, 
where G appears brightest, we determine that the center of the lensing galaxy is closer 
to the brighter s tellar component, a fact often seen among two-image lensed quasars (e.g.. 
Kayo et al.ll2010[ ). We are unable to determine a robust estimate of the Sersic index for G, 
due to the galaxy being blended with the quasar components. In order to extract magnitudes 
for G, we adopt the canonical value for elliptical galaxies, n = 4. This is justified by the fact 
that most galaxies acting as lenses are elli pticals, as spirals only con tribute to about 10% of 



the probability for multiple imaging (e.g., iFukugita fc Turner! 1 19911 ). 



1 



4. Laser Guide Star Adaptive Optics Imaging 



While the spectroscopic and imaging follow-up observations presented in the previous 
section indicate that SDSS J1334-I-3315 is a true lensed quasar system, its small angular 
separation makes it difficult to study the system in detail solely from these observations. 
One way to overcome this problem in the near-infrared is to obtain high-resol ution adaptive 



optics images, as already done for several other quasar lens sy stems (e.g., iMcKean et al. 
200?! : lAuger et al.lbood : Isiuse et al.lbood : Ihagattuta eraPboioh . 



We obtained high-resolution near-infr ared images for this qu asar system using the In 
frared Camera and Spectrograph (IRCS; iKobayashi et al.l l2000l ) at th e Subaru telescope 



along with the Laser Guid e Star Adaptive Optics system (LGS+A0188: lHayano et al.ll2008 



2010l : iMinowa et al.ll2010l ). A0188 uses a curvature sensor with 188 control elements and a 
188 element bimorph deformable mirror (BIM188), and operates at the Nasmyth focus of 
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the Subaru Telescope. IRCS was used in the low resolution mode with a pixel scale (y.'0527, 
producing a 54" field of view. The observations were performed in an engineering run on 
2011 February 18, with excellent seeing ~ 0'.'3 — O'.'S. We used an i? ~ 13.4 mag star located 
57" from the target as a tip-tilt guide star (see Figured]). We obtained a set of 60 s individual 
exposures with a total exposure time of 10 min for J-band, 4 min for if -band, and 12 min in 
K'-hand. The measured Strehl ratios were ~ 8% in J-band, ~ 12% in if -band, and ~ 12% 
in isT'-band. The full width at half maximum (FWHM) of the PSF was ~ 0"10 - 0"15 in 
all three bands. Therefore, although the PSF is not quite diffraction limited (~ 0.06" for 
Subaru Telescope in the i^'-band), it is better than what would be obtained with the Hubble 
Telescope at diffraction limit, in the equivalent band (~ 0.22"). 



The standard star P272-D ( iLeggett et al.ll2006[ ) was observed in the same night and at 



the same airmass with the target, without adaptive optics. Both the target an d standard 



star frames were reduced with the IRCS IMGRED package (JMinowa et al.ll2005l ). The zero- 
points which we estimate from aperture photometry are J = 24.39 ± 0.01 mag/s, H = 
24.53 ±0.02 mag/s, and K' = 23.78 ±0.02 mag/s, respectively, where the quoted errors refer 
only to the differences between aperture photometry of the standard star in different frames. 

We initially tried to build the PSF on an isolated star ~ 10" from the target, in the 
direction opposite to the tip-tilt star. However, this left significant residuals when subtracting 
the quasar components, most likely due to anisoplanatism introduced by the adaptive optics. 
We therefore built the PSF by fitting a two-component model, consisting of the Moffat 
and Gaussian analytical profiles, to the relatively isolated fainter quasar component. Our 
GALFIT modeling suggests a best-fit Moffat profile FWHM of ~ 6 - 8 pixels (~ 0"3 - 0"4), 
which is indicative of the seeing disk, and a Gaussian profile FWHM of ~ 2 — 3 pixels 
(~ 0"10 — 0"15), corresponding to the PSF core. This method significantly improves fitting 
and reduces the residuals after subtracting best-fit models. 

We show the adaptive optics images in Figure HI and summarize the results we obtained 
in Table 121 Thanks to the much improved angular resolution, the lensing galaxy is clearly vis- 
ible even before subtracting the stellar components. The spiral arm of the nearby galaxy Gl 
is also clearly seen. The relative astrometry derived from the LGS±A0188 image is broadly 
consistent with the MOIRCS result shown in Table [H but the attached errors are now much 
smaller. There are small discrepancies in the relative positions of the quasar images, larger 
than the error bars attached to the MOIRCS measurements, which must be due to an inaccu- 
rate original measurement of the position of image A (made difficult by the close proximity of 
the lens galaxy). The LGS±A0188 measurement corrects the MOIRCS-determined values, 
because the image A and the lens are more clearly separated. In addition, small (sub-pixel) 
discrepancies between the MOIRCS- and LGS±A0188-determined positions of galaxies Gl 



and G2 should be due to better detected morphological features in the LGS+A0188 imag- 
ing. The quasar image separation is found to be ^ = 0'.'833 ± 0.002. We also confirm that 
the center of the lensing galaxy is closer to the brighter quasar component. In addition, the 
high-resolution images enable us to make a more reliable measurement of the shape of each 
galaxy, such as the half-light radius and the elhpticity. Although the goodness-of-fit estimate 
X^ is almost insensitive to the choice of fiducial Sersic indices 1 (exponential profile) or 4 (de 
Vaucouleurs profile), less residuals appear in the K'-hand after subtracting the lens galaxy 
modeled by the Sersic index 4. This suggests an elliptical galaxy, in accordance with our 
original assumption. 

We must note that the photometry (especially that of the lensing galaxy) differs signif- 
icantly from the MOIRCS result. This is presumably because of the PSF and magnitude 
zero-point uncertainties, which are hard to characterize for AO images. On the other hand, 
relative photometry of the lens - quasar images should be more accurate in the case of the 
AO observations, due to the higher resolution. Also, objects (especially the galaxies) appear 
systematically fainter in the AO images, which may be due to uncorrected light scattered 
at large distances. While we leave for future work more comprehensive and sophisticated 
analysis of the PSF and photometric calibration, we confirmed that our astrometric results 
are rather robust against these uncertainties. 



5. Photometric redshifts 



In order to check the validity of our assumption that the strong absorption aX z = 0.557 
seen in the quasar spectra is caused by the lensing galaxy G, we estimate a photometric red- 
shift for this galaxy based on the magnitudes in the five observed bands (J, z, J, H, Kg/K'). 
We also estimate photometric redshifts for the two n earby galaxy companions. We emp loy 
the publicly available HyperZ (JBolzonella et al.ll2000l ) and EAzY (JBrammer et al.ll2008[ ) al- 
gorithms, both of which produce redshift estimates based on template-fitting methods. For 
both algorithms, we u se as templates the obs erved mean spectral energy distributions (SEDs) 
of local galaxies from IColeman et al.l (jl980l herea fter CWW), extrap o lated into ultraviolet 
and near-infrared with the evolutionary models of iBruzual fc CharlotI (119931 ). In the case of 
EAzY, we use a Bayesian prior on the apparent magnitudes, which is implemented in the 
code. The prior is added to help prevent the redshift probability distributions from having 
multiple peaks, as template colors can be degenerate with redshift. 

Since the best magnitude estimates for the lensing galaxy are significantly different 
between the MOIRCS JHKs- and the LGS+A0188 IRCS J HK' -hands, we use both of 
these estimates separately. On the other hand, for the two nearby galaxies, where the 
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aperture magnitudes are fairly similar, we use the values obtained with LGS+A0188. As 
photometric redshift estimates require absolute magnitudes, we attach error bars larger than 
the ones quoted in Tables [T] and [2], to allow for uncertainties in the magnitude zero-points. 

We conducted initial photometric redshift fits with HyperZ, and checked the results of 
the best-fitted spectral templates with EAzY. The resulting photometric redshift probability 
distributions computed with EAzY are shown in Figure [5l and summarized in Table [31 For 
the lensing galaxy, the best-fitted templates obtained with HyperZ are E (elliptical) and Sbc 
(spiral). The results are consistent with the probability curves calculated by EAzY. Although 
we cannot conclusively discriminate between the two templates based on the morphologi- 
cal information obtained even with LGS-I-A0188, both templates yield redshift probability 
distributions consistent with the redshift estimated from the absorption lines, within la con- 
fidence interval when using the LGS-|-A0188-estimated magnitudes, and at 1 — 2cr confidence 
interval when using MOIRCS-estimated magnitudes. We therefore conclude that the lensing 
galaxy is likely to be at 2; = 0.557, and is certainly at a redshift lower than that of the quasar 
components, consistent with the gravitational lensing hypothesis. 

For the bright nearby galaxy Gl, the probability curve for the Sbc template (matching 
the disk-like morphology) is in good agreement with z = 0.557, making it very likely that 
this galaxy is located at the same redshift as the lensing galaxy. It is also in good agreement 
with the previous redshift estimate of 0.50 ± 0.07, from the SDSS database. The G2 galaxy 
has a very poorly estimated redshift, due to its faintness, and also appears to be located at 
a redshift much larger than G and Gl. 



6. Mass Modeling 

We conduct gravitational lensing mass modeling to check whether the observed image 
configuration can be reproduced with reasonable mass models, as a final che ck of the gra vi- 



tational lensing hypothesis for this system. We use the public software glafic ( lOgurill2010l ) to 
solve the lens equation and to find best-fit mass models. We use 7 observational constraints 
for mass modeling: the (two-dimensional) positions of the lensed quasar images and the lens- 
ing galaxy, as well as the flux ratio between the quasar images. The positional constraints 
are taken from the accurate astrometric results derived from Subaru LGS+A0188 imaging 
(Table|2]). We use a flux ratio constraint of B/A = 0.75±0.05, consistent with most follow-up 
photometric results as well as the spectroscopic flux ratio over a wide wavelength range. The 
lens redshift is assumed to be z = 0.557. Given the small number of observational constraints 
available, we are forced to consider the simplest mass models: the Singular Isothermal Ellip- 
soid (SIE) and the Singular Isothermal Sphere with an external shear (SIS-I-7). Both mass 
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models have 7 parameters (the position of the lens and its velocity dispersion a or Einstein 
radius REin^ the ellipticity e or shear 7 and their associated position angles 6^, or 6^, and the 
position of the source), which indicates that there are no degree of freedom, as is common 
in the mass modeling of two-image quasar lenses. 

We find that both models can fit the lensed system perfectly, signifying that the choice 
of the models is reasonable. Table |4] summarizes the best-fit values for both models. We 
find that the best-fit ellipticity for the SIE model, e = 0.235, is significantly smaller than the 
observed elhpticity of the lensing galaxy, e ~ 0.62, but the position angles agree quite well 
with each other. Such large difference of the ellipticity and the agreement o f the position 



angle has often been seen in previous analyses of strong lens systems (e.g., iKeeton et al. 



19981 ). The position angle of the best-fit SIS+7 model is almost the same as that of the SIE 
model. In fact the angle corresponds to the position of the nearby galaxy Gl, suggesting 
that the quadrupole of the lens potential may be dominated by perturbation due to Gl. 
For both models, the total magnification and time delay are predicted to be /itot ~ 6.0 and 
At ~ 10 — 11 days, respectively. 

Here we discuss whether galaxy Gl can have a significant impact on the mass modeling. 
Based on the analysis in Section O we assume that the redshift of galaxy Gl is the same as 
that of the lensing galaxy. The follow-up photometry indicates that Gl is at least ~ 1.5 mag 
brighter than the lensing galaxy. Adopting the Faber- Jackson relation, we infer that the 
Einstein radius of Gl is at least twice as large as that of the lensing galaxy, i.e., i?Ein ~ O'.'S. 
Given the distance between G and Gl of d ~ 4'.'1, we estimate the external shear produced 
by Gl as 7 ~ i?Ein/2(i ~ 0.1, which is very close to the best-fit external shear of 7 ~ 0.096 
for the SIS-I-7 model, suggesting that Gl can have a significant impact on this lensed quasar 
system. A caveat is that Gl appears to be a disk-dominated spiral galaxy (see Figure [Hand 
Table [2]), and hence the use of the Faber- Jackson relation may not be appropriate. 

To further explore the origin of the quadrupole, we compute the goodness-of-fit estimate 
X^ as a function of the position angle. Specifically we compute x^ for each fixed value of 
the position angle, with the other model parameters varied to achieve a minimum x^. We 
find that the results are very similar between the position angle of the ellipticity ^e, for 
the SIE model, and the position angle of the external shear 9^, for the SIS+7 model, as a 
consequence of the well-known shear-ellipticity degeneracy. The result shown in Figure [6] 
indicates that the position angle inferred from mass modeling is perfectly consistent with an 
external shear from galaxy Gl, as well as consistent with the position angle of the elliptical 
lensing galaxy G at la level. On the other hand, the location of the nearby galaxy G2 
corresponds to an external shear with a large x^? which implies that G2 is not affecting 
the lens potential. While we cannot decompose contributions from G and Gl at this point. 
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detailed observations of a lensed host galaxy in deeper adaptive optics imaging observations 
may help clarify the origin of the quadrupole. 



7. Summary 

We have reported the discovery of the subarcsecond {9 = 0'.'833) gravitationally lensed 
quasar SDSS J1334+3315. The system has been identified as a lensed candidate from the 
SDSS data, using the standard SQLS candidate selection algorithm, despite the image sep- 
aration being smaller than the criterion for which the morphological selection is complete. 
Our follow-up observations at the Subaru and UH88 telescopes have confirmed the system 
to be a real strong lensing event. It consists of two images of a quasar at z = 2.426, lensed 
by a foreground galaxy. From the colors as well as absorption lines in the quasar spectra, 
we infer the redshift of the lensing galaxy to be 2; = 0.557. In particular, we have ob- 
tained high-resolution near-infrared images of this lensed quasar system using the Subaru 
LGS-I-A0188 system, with a better resolution in the i^'— band than would be provided by 
the Hubble space telescope in the equivalent band. The images clearly reveal the presence 
of the lensing galaxy in between the quasar images, and enable us to derive very accurate 
relative astrometry, as well as shapes of galaxies. In particular, the higher resolution made 
possible to identify errors in the non-adaptive optics astrometry, caused by the lens galaxy 
and one of the quasar images being blended, and thus made possible a more accurate mass 
modeling. We have found that the image configuration can well be reproduced with the 
standard mass models. We have pointed out that the nearby galaxy Gl, which is located 
< 5" from the lensing galaxy and is likely associated physically with the lensing galaxy, may 
affect the lens potential. This lensed system has the smallest image separation among the 
newly discovered SQLS lenses. 

This work represents the first observation and analysis of any gravitationally lensed 
quasar system with the Subaru LGS-I-A0188. Our results demonstrate the usefulness of laser 
guide star adaptive optics imaging observations, in particular that of the Subaru Telescope 
LGS+A0188 imaging, which has a much improved sky coverage, for the study of strong lens 
systems. As explicitly shown, the high-resolution imaging enabled by the adaptive optics 
system is crucial for accurate and robust measurements of the position and shape of the 
lensing galaxy, particularly for small-separation lensed quasars as presented in the paper. 
Indeed we are conducting an imaging survey of SQLS lenses with Subaru LGS-I-A0188, 
aiming to derive accurate astrometry and to detect lensed quasar host galaxies, which will 
further increase the value of the SQLS lensed quasar sample. 
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Table 1. Result of imaging observations with UH88/Tek2k and Subaru/MOIRCS 



Name 


AX [arcsec] 


AY [arcsec] 


I 


z 


J 


H 


Ks 


A 


0.000 ±0.003 


0.000 ±0.003 


19.16 ±0.03 


19.18 ±0.03 


18.57 ±0.02 


18.10 ±0.02 


17.41 ±0.01 


B 


-0.502 ±0.003 


0.673 ±0.003 


19.50 ±0.02 


19.75 ±0.03 


18.99 ±0.02 


18.40 ±0.01 


17.72 ±0.01 


G 


-0.159 ±0.015 


0.235 ±0.020 


21.50 ±0.37 


21.92 ±1.69 


20.29 ±0.11 


19.30 ±0.07 


18.53 ±0.06 


Gl 


-3.160 ±0.004 


-2.587 ±0.004 


19.95 ±0.04 


20.33 ±0.07 


18.68 ±0.03 


17.70 ±0.02 


17.06 ±0.02 


G2 


1.160 ±0.011 


2.073 ±0.013 






22.18 ±0.18 


20.42 ±0.08 


19.69 ±0.07 



Note. — The positive directions of AX and AY are defined by West and North, respectively. The astrometry is 
measured in the i^Tg-band image. The quoted errors include statistical errors only. The magnitude measurements for 
galaxies Gl and G2 arc photometric aperture magnitudes, whereas those of components A, B, and G arc model estimates 
fitted with GALFIT. The quoted errors include statistical errors only, and exclude uncertainties inherent in the PSF 
estimates or the zero-points. Magnitudes are given in the Vega system, without the Galactic extinction correction, o^ 

except for z-band magnitudes which are in the AB system 



Table 2. Result of Subaru LGS+A0188 imaging 



Name AX [arcsec] AY [arcsec] n Re [arcsec] e 6^ [deg] J H K' 

A 0.000 ±0.001 0.000 ±0.001 ••• •■• ••• ••• 18.71 ±0.01 18.27 ±0.01 17.53 ± 0.01 

B -0.484 ±0.001 0.677 ±0.001 •■• ••• ••• •■• 19.03 ± 0.01 18.54 ±0.02 17.83 ± 0.01 

G -0.134 ±0.006 0.235 ±0.007 =4 0.20 ± 0.03 0.62 ± 0.08 -42.7 ±5.9 21.01 ±0.06^ 20.06 ± 0.05^ 19.36 ± 0.05 

Gl -3.178 ±0.004 -2.578 ±0.002 1.68 ±0.06 0.50 ± 0.01 0.53 ± 0.01 -58.2 ± 0.7 18.68 ±0.04 17.90 ± 0.03 17.28 ± 0.03 

G2 1.159 ±0.005 2.033 ±0.004 2.46 ±1.09 0.08 ± 0.01 0.61 ±0.11 38.3 ± 7.0 22.69 ±0.31 20.52 ±0.13 19.94 ±0.10 



Note. — The positive directions of AX and AY arc defined by West and North, respectively, and the astrometry is estimated from the K'-hand 
image. Shapes of the galaxies (the Sersic index n, the half-light radius i?e, the ellipticity e, and the position angle 9e measured East of North) 
are all measured in the K'-h&nA image. The Sersic index of galaxy G is considered to be the fiducial value of 4, as explained in the text. The 
magnitude measurements for galaxies Gl and G2 are photometric aperture magnitudes, whereas those of components A, B, and G arc modcll 
estimates fitted with GALFIT. The quoted errors include statistical errors only, and exclude uncertainties inherent in the PSF estimates or the-j 
zero-points. Magnitudes are given in the Vega system, without the Galactic extinction correction. 

'^The quoted values of the lens galaxy magnitudes in the J- and H- bands correspond to the best-fit models (Sersic index 1 rather than 4). 
When constraining an index of 4, the values become ^0.1 mag lower.) 
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Table 3. Best-fit photometric redshifts 



Object & template 


Best-fit z 


\g limits 


2(7 limits 


G (AO), Sbc template 


0.694 


0.494- 


- 0.876 


0.329- 


-1.174 


G (AO), E template 


0.586 


0.397- 


- 0.768 


0.255- 


-0.959 


G (w/o AO), E template 


0.694 


0.530- 


-0.850 


0.368- 


- 1.004 


G (w/o AO), Sbc template 


0.818 


0.644 - 


-0.986 


0.456- 


- 1.202 


Gl, Sbc template 


0.513 


0.390 - 


- 0.642 


0.298 - 


-0.757 


G2, all CWW templates 


2.124 


1.774- 


- 2.579 


1.311- 


- 2.884 



Table 4. Best-fit mass models 



Model a [kms ^] e or 7 Qf, or B^ [deg] -REin [arcsec 



Mtot 



At [day] 



SIE 

SIS+7 



147 
143 



0.235 
0.096 



-47.1 
-47.2 



0.40 
0.38 



6.0 
6.1 



11.1 
10.0 
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Fig. 1. — Finding chart of SDSS J1334+3315 from the SDSS i-band image. The size is 
2'5 X 2'5. North is up and East is left. The bright star indicated by an open square was 
used as the tip-tilt star for the laser guide star adaptive optics imaging (see Sglfor details). 
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Fig. 2. — Follow-up spectra of SDSS J1334+3315 A and B taken with FOCAS at the Subaru 
Telescope. Quasar emission lines redshifted to Zs = 2.42 are indicated by vertical dotted 
lines. Strong absorption systems aX z = 0.557 and 2.172 seen in both component A and B 
are marked by vertical solid bars. The lower panel plots the ratio of the two spectra. 
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Fig. 3. — Follow-up images taken with Tek2k at the UH88 Telescope and MOIRCS at the 
Subaru Telescope. The size of each panel is 9" x 9". North is up and East is left. From 
top to bottom panels, we show /- and z-band images from Tek2k and J-, H-, and Kg- 
band images from MOIRCS. Left, middle, and right panels display original images, images 
after subtracting all but the lensing galaxy (component G), and images after subtracting all 
components, respectively. 
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Fig. 4.— Subaru LGS+A0188 images of SDSS J1334+3315 in J-, E-, and /sT'-bands. The 
size of each panel is 9" x 9". North is up and East is left. As in Figure [3l from left to 
right panels we show original images, images after subtracting all but the lensing galaxy 
(component G), and images after subtracting all components, respectively. 
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Fig. 5. — Photometric redshift probability distributions, normalized to unit area, for the 
three galaxies G, Gl and G2, calculated using the magnitudes in Tables [H and El In the case 
of G, we consider separately the IRCS JHK' (with Adaptive Optics) and MOIRCS JHKg 
(w/o AO) magnitudes. The probability distributions are calculated with EAzY, for the most 
consistent template estimates (see f|5]for details). The redshift of the quasar components, as 
well as that of the absorption line system measured in both quasars at z = 0.557, are also 
marked with vertical lines. 
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Fig. 6. — The best-fit x^ as a function of the position angle measured East of North. Since 
the results are very similar for SIE and SIS+7 models, the position angle can be interpreted 
as either Qf, of SIE or Q^ of SIS+7. Note that the other model parameters are marginalized 
over. The horizontal line indicates the \a range of the measured position angle of the lensing 
galaxy G. The vertical dotted lines show directions of the external shear corresponding to 
the location of nearby galaxies Gl and G2. 



